Effect of electron-phonon coupling on energy and density of states renormalizations of 

dynamically screened graphene 
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Electronic screening strongly renormalizes the linear bands which occur near the Dirac crossing 
in graphene. The single bare Dirac crossing is split into two individual Dirac-like points, which are 
separated in energy but still at zero momentum relative to the K-point. A diamond-like structure 
occurs in between as a result of the formation of plasmarons. In this work we explore the combined 
effect of electron-electron and electron-phonon coupling on the renormalized energy dispersion, the 
spectral function and on the electronic density of states. We find that distinct signatures of the 
plasmaron structure are observable in the density of states with the split Dirac points presenting 
themselves as minima with quadratic dependence on energy about such points. By examining the 
slopes of both the density of states and the renormalized dispersion near the Fermi level, we illustrate 
how one can separate fc-dependent and oj-dependent renormalizations and suggest how this might 
allow for the isolation of the renormalization due to the electron-phonon interaction from that of 
the electron-electron interaction. 

PACS numbers: 73.22.Pr, 79.60.-i, 73.40.Gk 



I. INTRODUCTION 

Graphene, a monolayer of carbon atoms, has been 
studied extensively since it was isolated in 2004»i Dis- 
coveries such as a giant Faraday rotation^ and large 
strain fields which are mathematically equivalent to elec- 
trons under 300T magnetic fields 3 - identify graphene as a 
playground for studying fundamental physics. Not least 
among these discoveries is the recent observation of plas- 
maronic peaks in the experimental angle-resolved photoe- 
mission spectroscopy (ARPES) data of Bostwick et ali^ 
for graphene epitaxially grown on H-SiC. Plasmarons^ 
are generally long lived quasiparticles arising from the 
coupling of the particle-hole continuum to charge density 
oscillations (plasmons)j 7 - They had not previously been 
observed in ARPES, but had been seen through other 
techniques in massive electron systems such as GaAs 
thins films and single crystal bismuthj&£ 

In the most recent ARPES work on graphene the sys- 
tem is doped by the addition of potassium giving rise to a 
diamond-shaped feature characterized by three energies 
labeled Eq, E\, and E^ (see Fig. [1]), which, for positive 
doping, we take as being implicitly negative as measured 
from the Fermi level. These structures were shown to be 
described by calculations of electron-electron interactions 
(EEI) in graphene via a dynamically screened^ random 
phase approximation (GoW-RPA)^^ Indeed the G W- 
RPA results show good agreement with experiment with 
the exception of a disagreement at high binding energies 
where Go W-RPA predicts two separate bands while the 
experiment shows these bands merging. A discussion of 
plasmarons in bilayer graphene is given in Ref. [l3|. 

Of course there are other possible interactions which 
could affect the electronic structure of graphene. Per- 



haps the most obvious is the electron-phonon interaction 
(EPI) which has been calculated from first principles, 
and is found to be dominated by optical phonons which 
occur near 200 meV.— In Fig. Q] we illustrate schemati- 
cally the primary modulations of the Dirac cone surface 
due to various renormalizations discussed in this paper, 
in order to introduce the characteristic energy scales and 
essential change in the Dirac cone structure. The renor- 
malization feature of a phonon at energy u>e and also the 
renormalization expected from a Go W-RPA calculation 
of the screened Coulomb interaction are shown. Here, 
rather than a simple renormalization which would mod- 
ify the effective Fermi velocity in the case of the EPI, 
the plasmaronic side bands (not fully shown) create a 
dominant diamond-shaped feature which signifies new 
band crossings. The inclusion of both interactions has 
a subtle influence on the observable plasmaronic struc- 
tures. The EPI results in (1 + A) renormalizations, in- 
cluding one causing a chemical potential shift, where A 
is the electron-phonon renormalization parameteri 15 i 16 
While for an ordinary metal, this (1 + A) factor modifies 
the effective mass, in graphene, the apparent result is a 
modification of the velocity of charge carriers, directly 
seen as a change in slope of the dispersion at the Fermi 
level. Electron-phonon effects provide other signatures 
beyond (1 + A) renormalization factors. For coupling to 
an Einstein phonon of energy loe there are structures in 
the density of states and kinks (illustrated in Fig. [1]) in 
the renormalized dispersion curves corresponding to this 
energy ! 15 : 16 There is also a background in the optical con- 
ductivity with sharp onset at we due to phonon- assisted 
absorption, separately revealed when twice the chemical 
potential is larger than uje ^ This is quite distinct from 
other possible effects such as band structure changes due 
to, for example, bilayers 1 -^, or from Landau quantization 



2 



under an external magnetic fieldJ^ In this last case the 
phonon peaks in the density of states becomes further 
modified.— 

While it is straightforward to identify features of in- 
dividual interactions it is often not clear, even to low- 
est order, what joint impact they have on the system. 
Experimental examination of the electronic structure of 
graphene provides information as to the overall renormal- 
izations but there is generally no method to disentangle 
the impact of different interactions. The issue becomes 
what 'knobs' exist in the system that can be controlled 
to distinguish these interactions. In the case of the EEI, 
the coupling strength (denoted by a, analogous to the 
fine structure constant of quantum electrodynamics) is 
directly dependent on the substrate dielectric material. 
Also, the self energy itself scales with doping. The EPI 
instead has a fixed energy scale (an Einstein frequency, 
oje) resulting in distinct features when the chemical po- 
tential is tuned above or below uje- Here, we show ad- 
ditional ways of extracting information on the EPI and 
EEI in graphene. 

In the following, we relegate to the appendix the pre- 
sentation of the details of the G W-RPA theoryi&±i>22r^ 
which sums the polarization bubble to obtain the EEI self 
energy. There we also present the self energy for an elec- 
tron interacting with a 200 meV phonon within a simple 
model.— 

In Section [Til we a Pply these self energies and show 
both the renormalized energy dispersion as well as the 
spectral function for a variety of EEI and EPI coupling 
strengths. In the absence of the phonon interaction, for 
strong a values the GoW-RPA theory shows a strongly 
displaced band as compared to the bare band at high 
binding energies (uj < E 2 ). At energies above the cross- 
ing at E\ (uj > Ei) the band shows much more modest 
changes. As a is reduced the strongly displaced band at 
high energies moves rapidly towards its non-interacting 
value. Further, we find that the dispersion and spectral 
function are modified substantially by the addition of the 
phonon self energy. In particular, the phonon coupling 
appears to restore weight to the bare band at high bind- 
ing energies due to the fact that its self energy has the 



where T, s (k, uj) is the self energy of cone s [we absorb the 
shift in chemical potential due to renormalization into 
the ReE s (fc,w)], ej? = svpk — fio is the bare dispersion 
and fio is the bare chemical potential. Note that we have 
taken h = 1. Also, we have removed the explicit vector 
notation from the momentum, fc, in Eqn |T]). In the cone 
approximation the dispersions, and therefore the spectral 
function, are not dependent upon the direction of k, only 



opposite sign to that of the EEI in this region which is 
more in agreement with what is seen in experiment. 

In Section IIII| we calculate the electronic density of 
states (DOS) with a focus on signatures of the plasma- 
ronic diamond structure in this quantity. We find that 
the Dirac-like crossings produce parabolic features simi- 
lar to those expected for a Dirac point when damping is 
included in the calculations! 1 ^ 16 These features shrink in 
amplitude for decreasing bare chemical potential, fi$, and 
are subject to broadening due to the EPI for |^ | > \^e\- 
Thus, there is a narrow region of doping in which one 
might hope to see features of the renormalized Dirac 
crossings in DOS based measurements. 

In Section Hvl we identify that the value of the density 
of states at the Fermi level is renormalized from its bare 
value due to the fc-dependence of the electron-electron 
self energy The renormalized dispersion and DOS have 
slopes near the Fermi level with dependence on both k 
and uj derivatives of the self energy. We show that there is 
a factor difference between these two slopes which allows 
for the separation of k and uj dependencies of these renor- 
malizations given one has information about both the 
DOS and renormalized dispersion near the Fermi level. 
It also allows us to separate, approximately, these two 
renormalization effects. It should be noted that, in this 
work, we restrict our discussion of renormalizations to the 
screened Coulomb interaction (which has both k and uj 
dependence) and the phonon interaction (which is taken 
to have only uj dependence). 

II. PLASMARON FEATURES IN A(k,uj) AND 
RENORMALIZED DISPERSIONS 

Angle-resolved photoemission spectroscopy (ARPES) 
allows one to probe the charge-carrier spectral density, 
A(k,uj), for a given momentum, fc, and energy, uj. The 
total spectral function in graphene can be defined as the 
sum of two components, one for the upper band, s = +1, 
and one for the lower band, s = — 1. We denote these 
contributions as A + and A~ respectively. This results in 
a total spectral function about a K-point 
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on its magnitude, k. 

We consider in this work two contributions to the self 
energy which we detail fully in the Appendix. There we 
introduce the electron-electron coupling strength, a, and 
the electron-phonon coupling strength, A, and take the 
total self energy for a given band, s, as the sum of these 
two interactions. Thus the total self energy for band s 
is £ s (fc,w) = T, EEI (k,uj) + >J EPI {uj,uj E ), where Y, EEI is 
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FIG. 1. (Color online) Schematic illustration of the renormalizations of the Dirac cone by the electron-phonon interaction, 
electron-electron interaction and the combined electron-phonon and electron-electron interaction. The EPI creates a feature 
at We- The screened EEI produces two crossing features at Eq and E2 with a plasmaron ring between at E\. Here, tJd is the 
energy of the Dirac point and cf is the Fermi level. 



the contribution of the EEI for the s band which is de- 
pendent on both k and w, and Y, EPI is the contribution 
from the EPI, which depends only on frequency u> and 

assumes an Einstein phonon mode at energy oje- The 

2 

electron-electron coupling strength is given by a = ^ v 
where g = g s g v is the combined spin-valley degeneracy 
factor, e is the electron charge, £0 is the effective dielec- 
tric constant and vp is the Fermi velocity. 

As has been noted in the literature,— for the case of the 
EEI there exists a diamond shaped feature in the renor- 
malized dispersion which has two band energy crossings, 
replacing the single Dirac point of the bare case, plus a 
plasmaron ring. Clearly, in the limit of a — > (no EEI), 
the GoW-RPA calculation must produce the bare conical 
dispersions. This is pertinent since the effective a goes as 
£0 , the inverse of the effective substrate dielectric con- 
stant, which is the average of that of the materials above 
and below the 2D graphene sheet. In this way, the EEI 
in graphene becomes tunable through the substrate. 

The renormalized lines are obtained from solutions to 
the equation 



ReE s (fc,w) 
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which corresponds to poles in the spectral density, 
A a (k,ui), in the limit of damping going to zero. We de- 
fine the renormalized dispersion as Et = e| + Rcl] s (fc, w) 
for a given band, s, such that ui — Ef, = as in Eqn. (|2|). 
In Fig. [2ja) we compare the bare and EEI renormalized 
bands for the case of a = 2. Here the solid blue (grey) 
curve applies to s = +1, the upper modified Dirac cone, 
and the red (light grey) curve to s = — 1, the lower one. 
Note that the axes are dimcnsionless and scale with dop- 
ing, here described through the chemical potential fi and 
the Fermi momentum, hp. 
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FIG. 2. (Color online) (a) The bare dispersion, et (solid black 
lines), and EEI renormalized energies E+ (blue) and E~ (red) 
for the case of a = 2. The upper-half cone is well approxi- 
mated by a simple renormalized cone (dotted line), (b) The 
total spectral density A(k,w) = A + (k,uj) + A~ (k, ui) with (c) 
A + (k,co) and (d) A~(k,u)) shown separately. 



For most momenta, Ef is reasonably well represented 
by a new cone (dotted line) with a single scaled renor- 
malization corresponding to a slightly steeper slope with 
the new Dirac point displaced to slightly lower energies, 
denoted by Eq. This is an energy regime of weak renor- 
malizations resulting in a dispersion which is close to the 
bare one. In addition, the fc-dependence of Y. EEI pro- 
vides multiple solutions to Eqn. ^ for a given u, the 
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FIG. 3. (Color online) Energy renormalizations and spectral 
functions for EEI with a fixed ot = 0.5 and varied EPI coupling 
A for /io = leV. (a-c) As the phonon strength increases, the 
renormalized band merges with the bare band (solid black) 
at high energy. The real part of the self energy due to EPI 
is large in this energy region and has opposite sign from that 
due to the EEI. Plasmaron structures are also modified, (d-f) 
The spectral function for the corresponding dispersions in the 
region of the plasmaron feature. 



results of which form a complicated bat-ears-like struc- 
ture in in the range —1.5 < cj/fio < —1-0. These have 
been identified in the literature as representative of plas- 
maronic effects^— The bottom of this structure ends in a 
point which matches the top of the red curve correspond- 
ing to the lower band, E^ , together forming a second 
Dirac crossing at k = and energy E2, below which the 
red curve again roughly has a conical shape, but this cone 
is far displaced from the bare cone for s = — 1 (solid black 
curve) an indication of strong Coulomb renormalizations. 



This band, E k , also provides an additional feature which 
peaks at Eq. 

Fig- Hlb) shows the spectral function and Fig^c) and 
(d) its sub-components. The full spectral function, shown 
in Fig. [2jb) , shows more clearly the appearance of new 
plasmaronic side peaks extending between oj = and 
E2, and then beyond but very broadened. Primarily one 
notes that the Dirac point has split along k = into 
two separate energies at Eq and E2. Comparing the 
renormalized dispersions of Fig. d^a) to the color plot 
of Fig. mb) we note the familiar diamond shape. This 
diamond is formed as a joint feature of the renormalized 
and plasmaronic bands. Below E2, the cone is strongly 
renormalized as this is an energy regime where transi- 
tions occur through coupling with the particle-hole con- 
tinuum and plasmons. In this energy range the spec- 
tral function again shows an additional peak not present 
in the renormalized energy band. However, comparison 
with the energy dispersions suggests that these peaks at 
larger k are the residual spectral weight along the simply- 
renormalized Dirac cone (dotted black). 

An interesting aspect of the color plots of Fig. [5] is 
to be found in the detailed examination of the variation 
in the spectral intensity as a function of k for fixed oj. 
For frequency just above Eq, for instance, it is clear in 
Fig. HUa) that there are three zeroes from Eqn. yet 
the intensity at these three values is quite different for 
each zero. This is a result of the /c-dependence of the 
imaginary part of the self energy. This is characteristic 
of the EEI and is in contrast to the EPI where the self 
energy has no k dependence, only an energy dependence. 
That the damping varies with k is clearly seen in Fig.[2jb) 
where we note a solid black 'x' shape, plus wings that are 
green, indicating an increase in smearing. 

We now consider the case of including the EPI. In 
Fig. [3] we display Ek dispersions and A(k, u>) color map 
plots for the a = 0.5 case for varied EPI coupling 
strengths, A = 0, 0.2 and 0.25 [(a)->(c)]. The black 
lines define the usual Dirac cone without interactions. 
In Fig. G^a), by comparing the solid blue curve with 
Coulomb correlations for a — 0.5 with the similar curve 
in Fig. [DJa) we see that reducing a has significantly 
changed the shape of the plasmaron structures of Ek in 
the energy region around the two Dirac crossings but 
more importantly, below this region the bare and renor- 
malized dispersions differ less in energy, in comparison 
to the a = 2 case, reflecting the reduced effect of the 
Coulomb interaction. In Fig.^b) and|3Ic), we show how 
the bands are modified by the additional inclusion of the 
200 meV phonon for increasing coupling strength. We 
see the usual ARPES 'kinks' at lo = ±uje representative 
of coupling to a boson32 We can see that the plasmaronic 
bands have changed shape rather significantly, but most 
noteworthy is the now nearly unrenormalized band in the 
lower-half cone. The electron-phonon interaction acts to 
merge these at larger binding energies. The renormaliza- 
tions and chemical potential change of the EPI counter- 
acts that of the EEI. For negative values of uj the real 
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FIG. 4. (Color online) Spectral function and self energies for 
a = 0.5 case as in Fig. [3] (a) Variation with k at uj = — 0.8/io 
for A + (k,uj) including only EEL The corresponding real and 
imaginary parts of the self energy are also shown. The value of 
uj — ReE — ek from Eqn. <|2j is plotted (black circles) and shows 
three zeros, while the spectral function has only two obvious 
peaks, (b) Variation with uj at fixed k = &f of the self energies 
for EEI (solid red), EPI (dashed purple) and EEI+EPI (solid 
blue) cases for the s — — 1 band as in Fig.[3{c) and (f) as well 
as the values of the left side of Eqn. @ with and without the 
EPI (blue and red circles respectively). 
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FIG. 5. (Color online) Imaginary part of the s = 1 band 
self energy for a — 0.5, A — 0.25 eV, and fio = leV as in 
Fig-EJc). The EEI is shown in red (light grey), EPI in dark 
blue (dark grey) and the sum is shown as the dashed curve. 
Inset: Enlargement of region around uj/fio = —1 with arrows 
indicating the locations of the three Dirac-like crossings at 
k = in Fig. Oc). 



part of the phonon self energy is positive while the EEI 
contribution is negative. Thus, one expects that, for suf- 
ficiently strong phonon coupling, the joint EPI and EEI 
will result in energy bands nearly identical to the bare 
case. 



This issue is further illustrated in Fig. EJd)-(f), where 
we display the spectral function for the cases shown in 
Fig.[2Ia)-(c) for a doping of /i = leV, and a = 0.5. The 
diamond shaped feature is still present in the spectral 
function. The height of peaks which do not appear in 
the renormalized dispersion will appear lighter or darker 
as Lj — Ek gets further or closer from zero, while not being 
a pole, and therefore not appearing in frames (a) through 
(c). 

We emphasize the fc-dependence of the EEI self energy 
in Fig. |4ja) where we show variations with k of some 
quantities for a given frequency uj = — 0.8^o which is in 
the region where the dispersion E£ has additional plas- 
maronic peaks [as shown in Fig. Gjja)]. We also include 
the quantity uj — RcS — the left hand side of Eqn (jJJ, 
shown as black circles. This quantity has three zeroes at 
this frequency which corresponds to three distinct lines 
of zeroes of Eqn. @ in the curves of Fig. Eta) on either 
side of k/kp = 0. The zero for smallest k, leads to the 
sharp peak in A + (k, uj) because the imaginary part of the 
self energy, shown as the solid black line in Fig. Eta) is 
also small. This is not the case for the other two zeroes, 
the last of which has been identified as the plasmaronic 
side band, here appearing as a rather disperse peak in the 
spectral function near k = 0.35/cp. The second zero oc- 
curs at a point of large damping and does not show up as 
a peak in A + (k, uj). The rapid variation of ImS^ 1 (k, uj) 
in this region is also the reason that the second broad 
peak in the spectral density (solid red curve) is shifted 
from the position in k of the third zero. The partial 
cancellation between EEI and EPI self energy effects is 
detailed in Fig.^Jb) where we see that in the region of in- 
terest the real parts carry opposite sign. In this case the 
self energy depends not only on uj but also on momen- 
tum. £_(& = kp,uj) is shown (red) for the EEI only case 
as well as the sum of EEI and EPI contributions (blue). 
We see that below uj = —1.5 the sum is rather small 
and this explains why the renormalized and unrcnormal- 
ized curves approach each other. Further, at these large 
binding energies there is only a single zero of Eqn. 
(solid red points) which shifts location with increasing 
electron-phonon coupling (solid blue points) towards pos- 
itive energies. While this is true for smaller values of a, 
larger electron-electron coupling may favor multiple ze- 
roes even at higher binding energy despite the presence 
of phonons. This result of the GoW-RPA calculation 
deviates from experimentally observed ARPES spectra 
across a wide range of a values^ The exact source of 
this deviation, or required additional interaction has yet 
to be understood. 

The imaginary part of the self energy and how it is 
changed with the introduction of an electron-phonon con- 
tribution can also provide additional important informa- 
tion on the plasmaron structure. In Fig.[5]we show results 
for — Im£ + (fc = 0,w) for varied uj along the k = line 
which defines the Dirac points. The figure applies to the 
case a = 0.5 and A = 0.25 eV which is shown in Fig.JHJc) 
and (f). The solid red (light grey) is the EEI contribution 
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while the solid blue (dark grey) is the EPI contribution. 
The sum of the two is represented by the dashed black 
curve. Note that the EPI provides no damping in the 
interval from —lje to uje relative to the Fermi energy at 
u> = 0. The physics behind this feature is that at zero 
temperature an electron of energy u (assume u > for 
dcfiniteness) cannot decay by the emission of a phonon 
if bj < lje- Thus, around oj — only electron-electron 
interactions contribute to the damping. Another impor- 
tant feature of the EPI is that it gives, on its own, a zero 
in — ImE + (fc = 0, uS) at oj = —lue — Mo which occurs at 
— 1.2 eV in the case shown. The imaginary part is pro- 
portional to the density of states which has been shifted 
away from the Fermi level by ±uje [see Eqn. IA. 15] and 
therefore the zero in the bare quasiparticle density of 
states at — Mo appears as a zero in ImE at —lje — Mo- In 
this region of energy, the EPI adds little to the damp- 
ing which is mainly provided by the EEI which shows a 
strong peak just below u = — Mo- This peak damps out 
structure in this region of the energy dispersion as seen 
in Fig. [3Jc) where there appears to be three Dirac-like 
crossings at k = but the middle one does not appear as 
a Dirac point in Fig. [3Jf ) . This is further illustrated by 
the inset of Fig. [5] which locates with arrows the positions 
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where N = 

We present the results of Eqn. fl]) for the density of 
states using the GoW-RPA in the absence of the EPI in 
Fig-IH for varied a. The bare DOS (green dashed line) is 
a wedge shape pointing to the chemical potential, with a 
value of zero at mo in the absence of scattering. For non- 
zero a, this point is lifted due to the finite lifetime at this 
energy provided by | ImE |J£ Additionally, the slope of the 
DOS is modified; reduced in magnitude with increasing 
a over a broad energy scale. Furthermore, in Fig. |6](b) 
we examine the region near the Dirac point and observe 
distinct features in the range u> £ [ — 1.6//0, — Mo]- I n the 
DOS, the Dirac crossings of Ej, at E and E 2 result in a 
parabolic minima located at those energies as well as a 
similar feature at E\ associated with the plasmaron ring. 
It is these minima in the DOS that we use in this pa- 
per to identify the energies Eq, E\ and E 2 rather than 
the corresponding features in the renormalizcd dispersion 
curves. Plasmaron structures, while small in the DOS are 
in principle imprinted on N(u>) between the two minima 
associated with the split Dirac points at k = 0. For vari- 
ation in a we observe that these energies are modified. 
This is plotted in the inset of Fig. Mjo). The energies scale 
with /iq for a given a. The energies, E$, E\, and E 2 and 



of the three Dirac-like crossings from Fig.[3](c) and shows 
that the middle one lies at the energy of large damping 
as compared to the other two arrows which occur in a 
region of low damping. 



III. PLASMARON FEATURES IN THE 
DENSITY OF STATES 

The electronic density of states at a frequency u> is 
given in terms of the charge carrier spectral function by 

N(u)=g 8 g v J2M^u)> (3) 
k 

where g s and g v are the spin and valley degeneracies, re- 
spectively, and A(k, uS) is given by Eqn. ([1]). While the 
EPI is strongly dependent on the value of moj m the ab- 
sence of the EPI (A = 0), the resulting density of states 
including EEI is nearly independent of the value of the 
bare chemical potential with only a very weak depen- 
dence at large uj approaching the band cutoff, W c . In 
general we can write the density of states as 



(4) 
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their variation with a are of importance for the experi- 
mental community, as different substrates can modify a 
through several orders of magnitude, from 2.0 for H-SiC 
to 10~ 3 for SrTi03. In the latter case, there is evidence 
for strong variation of the substrate dielectric constant 
with temperature, creating a scenario where the value of 
alpha is modified by an order of magnitude as tempera- 
ture is changed from room temperature to liquid helium 
temperature^ 

In Fig. [7J wc display for the a = 2 case the DOS with 
(blue) [dark grey] and without (red) [light grey] a 200 meV 
phonon for chemical potential above (a) and below (b) 
uie- For too large a chemical potential, the DOS be- 
comes smeared in the region of E , Ex and E 2 due to 
additional broadening provided by the |ImE BP/ | result- 
ing from electron-phonon scattering. This is similar to 
the lifting of the Dirac point in the EPI only case for 
Mo > uje^ Decreasing the chemical potential to avoid 
this smearing, however, will weaken the plasmaronic fea- 
tures in the DOS due to the scaling with /io- As a result, 
there is perhaps a region of doping values just less than 
the phonon frequency where one might look for plasma- 
ronic features in tunneling. 

Several features in Fig. [7J are to be noted. The intro- 
duction of coupling to a phonon has shifted the structures 
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in N(lu) corresponding to Coulomb rcnormalizations (Eq, 
Ex, E 2 ) towards positive energy though the distance be- 
tween the three energies remains fairly constant. There 
are also additional structures introduced which corre- 
spond to pure phonon peaks seen at uj = ±lue = ±0.2 
cV. We note also that while the slope of the DOS curve 
at the Fermi energy, uj — 0, has changed, the value of 
the DOS at this point is unchanged by the introduction 
of the EPI. However, its EEI renormalizcd value is quite 
different from the bare density of states (dashed line) . 

Taking a derivative with respect to oj acts to accen- 
tuate features of the parabolic minima and has been 
shown in the past to be an effective way to bring out 
subtle structures in general. 15 ' 26 In Fig. [8ja), we show 
N'(uj) = -j^- for three values of /i for the EEI with 

a = 2. Also shown are the derivatives of the bare DOS, 
shown with dashed lines, which are smeared step func- 
tions at iM). In the EEI case, one immediately notes the 
suppressed value of N'(u) at u = 0. Further there are 
new features below —fio- These features include several 
zeroes and shrink in magnitude with reduced /Uo but do 
not change their relative shape which reflects the scaling 
of the self energy with /i noted for the pure EEI case. 
Seen most clearly for fig = 0.25 eV (black), there are five 
zeroes below —fAo- The zeroes which cross from negative 
to positive (having positive second derivative) correspond 
to the local minima in the DOS which are caused by the 
features in the spectral function at energies E 2 , E\, and 
Eq. Thus, these zeroes describe these special energies 
and represents a clear image in the quasiparticlc density 
of states of Coulomb renormalizations and particularly 
of plasmarons in graphene. If one includes the EPI as in 
Fig. [5th), there are now two new features, one occurring 
at lj = —ijJe and another at —fio - wb- This latter energy 
is where the llmE^ 7 ! goes to zero. As a result of the 
EPI the value of N'(oj) is brought closer to the bare case 
both at uj = and especially so at large negative u> due 
to the feature at — ^ — ^e- While in the EEI case the 
energies of the zeroes of these curves scale with /j, , in the 
EPI case this is not so. 



IV. RENORMALIZATIONS IN E k AND THE 
DOS 

The observed impact of the EEI and EPI on slopes of 
the dispersion, Ek, and the DOS at the Fermi level as 
well as their variation with coupling parameters, moti- 
vates the detailed study of these renormalizations. To do 
this we start by expanding the self energy about a given 
momenta, k = k* and energy, oj = uj* to lowest order. 
We write the real part as 

ReE(fc,w) ReS(fc* ,oj*)+-?-ReE(k,uj*)\ k ,(k - k*) 

ok 
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FIG. 6. (Color online) DOS for EEI only, (a) Increasing a 
results in reduced slope far from the Dirac point, (b) Closeup 
of upper frame. Increasing a results in distinct parabolic fea- 
tures with the minima positioned at the energies correspond- 
ing to the features seen in Fig. [2£b) which are denoted as -Bo, 
Ei and Ei- Inset: Eq, Ei and E2 versus a. 



We introduce the notation 
1 d 



\k=k' = —-^rRcT,(k,uj*)\ k , 
vp ok 

d 

^lo=w = -T^Re£(£;*,u;)|^. 
ou 



(6) 
(7) 



(5) 



From these we define k and oj dependent renormalization 
factors as Z k = 1 + Afe and Z u = 1 + A w . At the Fermi 
level, Eqn. now becomes uZ w — t^Z k = 0. From this 
we can approximate the renormalizcd energy in terms of 
the bare dispersion and the renormalization factors as 
E k = e fe f^. 

We show details for two points on the dispersion 
curves, the Fermi surface (k* = kp, oj* = 0) and the 
lower Dirac crossing (k* = 0, lj* = E 2 ). The slopes 
of Ek at these two points are indicated in Fig. [9jc) [for 
q = 2] as dashed red and dashed blue, respectively. The 
slope of Ek at the Fermi surface defines quite well the 
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FIG. 7. (Color online) DOS with EEI + EPI for a = 2, 
A = 0.25: (a)/io > oje, (b) fio < u>e. Phonon interaction 
shifts plasmaronic and Dirac features towards positive energy 
due to the shift in chemical potential. Inclusion of additional 
phonon scattering acts to fill between these features making 
them less apparent in the DOS. Insets: larger field of view of 
DOS. For reference we include the bare DOS (dashed line). 




-0.6 -0.4 -0.2 -0.6 -0.4 -0.2 
CO (eV) CO (eV) 



FIG. 8. (Color online) N'(u) = j^^ 1 for a = 2, A = 0.25 
with: (a) EEI and (b) EEI +EPI° for"/i = 0.25 (blue), 0.15 
(red), and 0.1 eV (black). Also shown are the bare cases 
(dashed lines). 
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FIG. 9. (Color online) Renormalizations for EEI in GoW- 
RPA for a range of a values shown at (a) the Fermi level and 
(b) E2- (c) Energy dispersion for a = 2 (solid line) com- 
pared to the simple renormalizations at u> = (red dashed) 
and E2 (blue dashed). The u — renormalization, shown 
schematically by the line with rise and run given by Zk and 
Z^ respectively, roughly agrees with the full calculation until 
nearly E\. 



variation of the Coulomb rcnormalizcd dispersion curves 
over the entire range down to the first Dirac point at 
E n . The dashed blue line which defines the slope of the 
renormalized dispersions at the second Dirac point, E2, 
begins to deviate somewhat from the solid blue curve for 
energies away from E2. Thus the extent over which this 
approximation scheme is valid is limited near the vicinity 
of E2 ■ It is worth noting that in general, the slope of this 
blue dashed line is not the same as that of the red dashed 
line. 

Fig. [3Ja) and (b) show slopes at the Fermi energy, 
w = 0, and E 2 , respectively, as a function of a. The 
solid (black) curve gives Z k , the dashed-dotted (red) , Z^ 
and the dashed (blue) Ek/et- Wc define this ratio to be 
Ek/tk = itt-c = Zk/Z U i which is rather close to one for 
all values of a even though for a = 2, Zk=k F ~ 1.7 and 
Zu)=q ~ 1.5. At to = the Zu and Z^ renormalizations 
both grow in tandem for increasing a. This is not the 
case for other energies, such as E2, where above some 
coupling strength the renormalization values cross. This 
is a consequence of the fact that the EEI energy renor- 
malization limits to zero at the Fermi energy, not at £'2- 

To illustrate this point, we define a quantity /3 = x eei 

k = k p 

plotted as the dotted black line of Fig. HJa). It is clear 
that /? is essentially constant for all values of a and in 
fact has a value of roughly 2/3. If one defined a similar 
quantity at E2 , or any other energy, this behavior would 
not be observed. In this light, /3 has a universal value for 
varied a, or to put it in a more experimentally relevant 
context, /? is substrate independent. 

This is an important result which needs to be kept 
in mind when analyzing ARPES data with a view at 
understanding the role of Coulomb interactions. Also 
shown, for completeness, is the chemical potential shift, 



9 



Afx/fxo, due to the EEL Note that the Zs begin at a 
value of 1 for a = 0. The lower part of the left frame gives 
similar results for the Dirac point at u = E^ . For the well 
known case of the electron-phonon interaction Zk = 1 
and we recover the known result that E). = efc/(l + A^ PI ). 

With this idea of slope renormalizations near the Fermi 
level, we draw the reader's attention back to the DOS 
calculations of Fig. [7] We can understand both the DOS 
value and slope change at the Fermi level. One must 

rewrite the delta function o(eZj. — ujZ u ) = — \z k \ k • 
Thus the density of states is subject to this Z^ renor- 
malization. 

Further, in the limit u — > 




— C E ""- (7) 

• • Eqn. (ll)for p=2/3 



0.5 



1.5 



N(u -» 0) _ N bare {co§^) N bare {0) 
N~ ~ \Z~Z\ \Z kF \ 



(8) 

where N bare (uj) = \u + (jlq\. We conclude from Eqn. flSJ) 
that the density of states at the Fermi surface for an ex- 
ample case of a = 2 where Z^ F = 1.63 provides about a 
40% reduction to the DOS as seen in Fig. E^a) at u> = 0. 
This is to be contrasted to the well known effect in metal 
physics that the electron-phonon interaction leaves the 
DOS at co = unaltered (Zf F P/ = 1). Further, while the 
variation in N(ui) for small w remains linear in lo, its slope 
is changed by the factor = Z u)= q/Z'^ f which is to 
be contrasted with our previous result for the dispersion 
curve renormalization where the renormalized quasipar- 
ticle energy, E\ = z kF ^ g = m e efc. We now have slope 
renormalization factors for the DOS and dispersion given 
by tun and m e , respectively. Measuring in scanning 
tunneling spectroscopy (STS) and m e in ARPES we can 
separate the renormalization factor, Zk F = (m e mjv) _1 - 
If we include both the EEI and the EPI we have that 



z k =i+xr T , 

z u = i + x EEI 



+ A, 



EPI 



(9) 
(10) 



FIG. 10. (Color online) DOS and Ek F slope renormalizations, 
mjj and m e at the Fermi level due to the EPI and EEI. Manip- 
ulation of these slopes using Eqn. (|11[) allows the extraction 
of the EPI contribution A5=o given knowledge of the factor 
/3 in the EEI. The result of Eqn. © and Eqn. (JTTJ match 
precisely due to the independence of the EEI on /iq. 



gle curve for the density of states whatever the value of 
chemical potential. As a consequence, the contribution 
to the renormalization at the Fermi surface from the EEI 
is constant for variation in /io- On the other hand, the 
electron-phonon mass renormalization increases with in- 
creasing /xo' 16 ' 17 As a result, when both the EEI and EPI 
are considered, m t and tun will vary as seen in Fig. 1101 
While in this theoretical work, the evaluation of Eqn. pT|) 
is trivial, we include the case for j3 = 2/3 to emphasize 
that the knowledge of m e and mjv is sufficient to ob- 
tain the EPI renormalization, A^£o, shown in Fig. [TO] 
This analysis should be an experimentally realizable task, 
given one performs ARPES and STM on a sample with 
the same a. 



If we note that (mjvm^ 1 = Z^a = 1 + X^Eo + 
we can obtain both Xff 1 and the sum of A^fg and A^o 
separately. Now the value found for Aj^f 1 can be used to 
deduce [from Fig.[9ja)] the value of the coupling constant 
involved (a) and consequently, the substrate dielectric 
constant. Also from examination of Fig.^a) we see that 
for a given a we can relate Af^f 7 to X^Eq through j3 such 
f3XE EI ■ To continue this example, we would 



that X^Eq 
find that (mjvw 
conclude that 



ft* 

2\-l 



/3X EEI + X e Eq horn which we 



X 



EPI 




P(m N m e )~ 



(1-/5). (11) 



This analysis should allow a direct experimental estimate 
of the strength of the electron-phonon interaction given 
m e measured by ARPES and rriN by STS. 

As we have seen, when only Coulomb effects are con- 
sidered in the GqW-RPA approximation there is a sin- 



V. CONCLUSIONS 

We have studied how Coulomb interactions change 
the density of states in graphene within the dynamically 
screened GoW-RPA. As was previously found for the dis- 
persion curves, a single universal function scaled by the 
chemical potential N \uj / no) / fio can describe all doping 
values for a given value of the coupling strength a. The 
Dirac point of the bare case splits into two as predicted 
and seen in ARPES experiments and these show up in 
the DOS as two minima at Eq and E2 with a quadratic 
behavior. The plasmaronic ring presents itself as a fur- 
ther additional minimum at energy E\ between Eq and 
E 2 . These structures can serve to identify plasmarons in 
the DOS as measured in STS experiments. Application 
of a first derivative to N(u>) further enhances these sig- 
natures. We describe how the three energies vary with 
coupling a which is inversely proportional to the sub- 
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stratc dielectric material. 

We find that Coulomb coupling reduces the value of 
the interacting DOS at the Fermi surface below its bare 
value by a factor of Zk F , independent of doping. This 
is a radical departure from the well known result that 
the electron phonon interaction itself does not renormal- 
ize the DOS at uj = 0. This difference is understood as 
a result of the fact that Coulomb interactions provide a 
self energy that depends on both momentum, k, and en- 
ergy, u, while the EPI depends only on u. The factor 
Zfe F comes directly from the derivative with respect to k 
of the self energy. On the other hand, the slope of the 
renormalized dispersion curves at Uf is given by a ratio 
Zk F /Z u= o, where Z^ is related to the derivative of the 
self energy with respect to uj. Both factors have simi- 
lar magnitudes which results in the ratio ranging from a 
value of 1 at a = to 1.2 at a = 2. 

In a similar fashion, we show that the slope of the 
linear-in-w dependence of the DOS at the Fermi surface is 
given by the factor Z u =o / Z 2 p . Consequently a joint mea- 
surement of the slope of the dispersion curves by ARPES 
and that of the DOS by STS allows one to measure di- 
rectly the Coulomb renormalization factor Z^ F as the 
ratio of these two slopes. 

We have also studied the effect of additionally includ- 
ing the electron-phonon interaction. This second inter- 
action can significantly affect the plasmaron structure 
in both renormalized dispersions and density of states. 
In particular, the EPI causes the energies Eq and E2 
to move towards the Fermi level as does E\ which re- 
mains roughly at the same relative position between the 
quadratic minima associated with the split Dirac points. 
For the set of parameters examined, an important effect 
of including both EPI and EEI is that well below the sec- 



Here the non-interacting Green's function in the s' band 
is given by 

G° s ,{k + q,oj + in m )= 1 — , (A.2) 

UJ + iQ. m - e% +q 

where i£l m = i2irm,m = 0,±1,±2... are the bosonic 
Matsubara frequencies and e k , is the energy, relative 
to the Fermi level, of the s' band at final momentum 
k' = k + q. Due to the symmetry of the cone approx- 
imation, the k direction can define the coordinate sys- 
tem, such that e k , is a function of the magnitudes of k 
and q and the angle between them, 9k q - This results in 
e l' = vfs' ^Jk 2 + q 2 + 2kq cos 6kq - /lo- 
in addition to an effective potential, the GoW-RPA 
calculation for a chiral system includes a scattering am- 



ond Dirac point, the real part of the self energy of each 
interaction separately carries the opposite sign, and so 
partially cancel against each other. This leads to a re- 
duction of the difference between renormalized and bare 
dispersion in the high binding energy region than would 
arise from the EEI alone. 

We also establish a procedure whereby the EEI and 
EPI rcnormalizations at the Fermi surface can be sepa- 
rately determined from combined measurements of the 
DOS in STS and dispersion curves in ARPES. 
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Appendix: Theoretical Background 

1. Electron-Electron Interaction (EEI) 

One can calculate the self energy within a random 
phase approximation (RPA) through a procedure which 
is standard for a 2-dimensional electron gas (2DEG) 20 
and modified to include the issue of chirality in a chi- 
ral 2-dimensional electron gas (C2DEG). The electron- 
electron interaction is assumed to be due to an effective 
potential, W, which accounts for a screened Coulomb 
repulsion. The form of the self energy then includes a 
Green's function, G, and the effective screened potential, 
W, leading to the GoW-RPA expression in terms of a 
sum over Matsubara frequencies, given by 



(A.l) 



I 

plitudc which is the overlap of the state k in the s band 
with the state k! in the s' band. This factor acts to re- 
move backscattering, and is given by 

F ss <(9 hk ,) = i[l + cos(0 fcfc O** / ], ( A -3) 

where Of-y is the angle between the vectors k and k' . 
Again, if one defines k' relative to the direction of k, we 
can write the angle between as 

flfe^arctanf -fj 11 ^ Y (A.4) 
\yk + q cos 9 kq J 

where q and k are dimensionless quantities as defined 
below. We note that in calculations, it is important to 
correct for the non-principle value of the arctan function. 



I 

Ef PA (fc, W ) = -T J2 G°,(k + q,u; + in m )F ss ,(6 kk ,)W(q,iSl m ). 
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The final piece to Eqn (|A.1[) is the effective potential 
given by 



W{q,iil m ) 



v a 



i-v q u ( q ,tn m y 



(A.5) 



where V q is the 2D Coulomb potential, V q = , where 
£o is the effective dielectric constant of the medium and 



IIo is the polarization function for doped graphenei 10 ' 22 
All necessary quantities can be written in a dimen- 
sionless form by scaling the energies by the bare chemi- 
cal potential, fj,Qi and scaling the momenta by the bare 
Fermi momentum. We denote these dimensionless quan- 
tities with a bar notation. With this in mind we define 



V q U 



a, where a 



[g = g s g v is the combined 



spin and valley degeneracy factor] and II is given for a 
given q and ^^2 as^ 2 - 



n(9,ifl) = -i- J 



8 V^ 2 + q 2 



=Re 



arcsm 



2 + ifl 



2 + iQ\ I 



(A.6) 



Alternate polarization functions have been explored such 
as, for example, the case of gapped graphene 2 ^ which 
maps onto Eqn. (|A.6|) as the gap goes to zero. In this 
case we can write the effective potential as 



W(q,in) 



2ira vf 



1 



g kFq — oW{q 1 in) 



(A.7) 



Evaluation of the self energy of Eqn. (|A.1|) is generally 
performed in two parts, identified as the line and residue 
(RES) portions^ The line component is completely real 
while the residue portion has both real and imaginary 
parts. These are given by 



OO 27T 



s'=±l 



dqd6 kq a _ x 



o o 



2tt g 



£ 1 (q,u - e s k+q )F ss ,(6 kk >)x 



(A. 



and 



oo 2tt 



s'=±l 







d qd9 kq a 
~2^T~g 



F ss '(O kk .')x 



2^ £ 



ui — e 



k+q 



(A.9) 



where £ 1 is given by 



£ = 



1 - V q Il q-aU 



(A.10) 



Here we have also written the band energies for general 
k and q in a dimensionless form 



= s'Jk 2 + q 2 + 2kqcos0 k „ -I (A.ll) 
Mo v 



which is the energy in the Dirac cone approximation. 



The integral over Cl in Eqn. (|A.9[) contains two terms, 
the first has a Lorentzian shape which provides a natural 
cutoff while the second term is odd and vanishes over the 
H integral. Also, n for purely imaginary frequencies is 
completely real.— The evaluation of real frequencies in, 
for example, £ _1 (g, u>) requires that one include a finite 
lifetime for the plasmons. This can be done in the form 
£ _1 (g, ui) — > £~ 1 (q, u> + iT), where T is the plasmon scat- 
tering rate 



10.28 



Failure to include such a term results in 
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infinitely long lived plasmons which cannot scatter with 
electrons. Thus, such a rate controls electron-plasmon 
coupling (plasmarons). 

Previous worki£ has shown that the g-integrals of 
Eqn. (|A.8|) and (|A.9|) suffer from a divergence at large q. 
This ultraviolet divergence is an artifact of the assumed 
linear dispersion in graphene, which allows for arbitrarily 
large q scattering. A detailed calculation would include 
a band cutoff at an energy W c ~ 7eV which corresponds 
to the energy which maintains the number of states in 
the Brillouin zone. Hence, this is important for density 
of states calculations, such as Eqn. where we include 
a normalized ultraviolet cutoff of A = ±W c /[J>o- A simi- 
lar cutoff is required in the self energy calculation itself. 
It has been established^ that the value of A lies in the 
range of 10 — > 100 for experimentally relevant values of 
chemical potential. Thus we have chosen a fixed cutoff 
of A = 50 in the self energy formulae, independent of the 
value of fiQ. This fixed cutoff removes any fj,o dependence 
in the EEI self energy. If one seeks to compare directly 
to experiment, then A can be modified for a given dop- 



for the real part, 



and is zero outside these intervals. We see that this 
self energy is a function of frequency, u, and also de- 
pendent upon the choice of Einstein frequency, u>e- Wc 
assume a model where = 200 meV as has been done 
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